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New exact results are obtained for relativistic acceleration of test positive ions in the laminar 
zone of a planar electron sheath evolving from an initially mono-energetic electron distribution. 
The electron dynamics is analyzed against the background of motionless foil ions. The limiting 
gamma-factor 7 poc of accelerated ions is shown to be determined primarily by the values of the ion- 
electron charge-over- mass ratio fi — m e Zp/m p and the initial gamma- factor 70 of the accelerated 
electrons. For fi > | a test ion always overtakes the electron front and attains 7 p0 o > 70- For 
(Jt < 5 a test ion can catch up with the electron front only when 70 is above a certain critical 



value 7 cr , which for /i< 1 can most often be evaluated as 7 C 



i/xexp(^i 1 



In reality the 



protons and heavier test ions, for which 7 cr > 10 is enormous, always lag behind the front edge 
of the electron sheath and have 7 poo < 70; for their maximum energy an appropriate intermediate 
asymptotic formula is derived. The domain of applicability of the laminar-zone results is analyzed 
in detail. 

PACS numbers: 52.30.Ex, 52.38.Kd, 52.40.Kh 



I. INTRODUCTION 

One of the most impressive latest achievements in 
laser-plasma interaction has been the observation of well- 
collimated high-energy proton and ion beams produced 
from thin metallic foils irradiated by ultra-intense sub- 
picosecond laser pulses pj 0, Q. An exceptional 
quality, demonstrated recently for thus produced proton 
beams |5j , makes them very promising for many potential 
applications 0. 

In its gross features, the mechanism of ion acceleration 
in the cited and other similar experiments is believed to 
be reasonably well understood, and has been nicknamed 
TNS A (target normal sheath acceleration) 0, ■ High 
directionality and low phase volume of generated protons 
indicate that they are accelerated in a highly ordered 
electric field normal to the virtually unperturbed planar 
rear surface of the laser-irradiated foil. The electric field 
is caused by charge separation in the sheath layer that is 
formed by energetic electrons generated by absorption of 
the laser pulse. 

However, an attempt to provide a more detailed and 
comprehensive theoretical description leads to a very 
complex system of plasma dynamics equations that can 
hardly be ever solved rigorously. To establish practically 
useful dependences and relations, one has to introduce 
additional simplifications. As is typical in other areas of 
physics, of particular value for gaining a deeper under- 
standing of the process of ion acceleration prove to be cer- 
tain particular idealized but exactly solvable problems. 
Salient and well known examples are (i) a self-similar evo- 
lution of the ion distribution function by quasi-neutral 
plasma expansion into vacuum and (ii) a virtually 
exact two-fluid solution of the isothermal plasma expan- 
sion with a full account of charge separation effects 
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In this paper we present a rigorous solution and a full 
parametric analysis of another such idealized problem, 
namely, the problem of test ion acceleration in a dynamic 
sheath of relativistic electrons with the delta- function ini- 
tial velocity distribution. 

Typically, fast protons in laser irradiated metallic foils 
originate from a thin (few nanometers) contaminant layer 
of water and hydrocarbons at the foil surface 0, 0, □J. 
Then, a natural simplification would be to assume that 
the heavy bulk ions (like Au for example) of a metal- 
lic foil are infinitely heavy and stay at rest, while the 
accelerated protons are treated as test positive charges 
initially located at the foil surface. Our present solution 
is essentially based on this assumption. 

Next, one has to choose how to treat the electrons. 
A widely used assumption is that the electrons instan- 
taneously relax to the equilibrium Boltzmann distri- 
bution in the time-dependent electrostatic potential of 
the expanding plasma: it was employed in both solu- 
tions 0, cited above. With immobile ions, such 
an assumption allows straightforward calculation of the 
electrostatic sheath potential <fi(x) either in the one- 
temperature |12( or multi-temperature [13j cases. An 
obvious problem with this approximation is that it leads 
to a diverging result for the maximum energy of ac- 
celerated test ions because the corresponding potential 
4>{x) = -(2T/e)ln[l + x/(-\/2A D )] El logarithmically 
diverges at x — ► +00 in the planar geometry; here Ajj 
is the Debye length at the base of the electron sheath 
with temperature T, +e is the elementary charge. How- 
ever, as was pointed out by Gurevich et ai. this 
divergence is not physical because even if one assumes 
that the hot electrons have a perfect Maxwellian distri- 
bution initially, at t = 0, it still takes an increasingly 
long time At ~ t eq (x) oc x for the Boltzmann relation 
n e oc exp(e0/T) to establish at an increasingly large dis- 
tance x from the initial plasma surface (for more details 
see section lTlI B 2l below') . As possible remedies, attempts 
have been made to use ad hoc quasi-equilibrium electron 
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distributions truncated either at high velocities 0, 0] 
or at large distances ^(| . Evidently, neither of these two 
approaches is fully self-consistent. 

Without the Boltzmann relation, a self-consistent 
treatment requires that one starts with a given initial 
electron distribution function at t = 0, and then calcu- 
lates its evolution for t > 0. For high-energy (multi-MeV) 
electrons this can be done in the collisionless approxi- 
mation. In this work we solve this problem in the sim- 
plest case of initially monoenergetic electrons, i.e. when 
at t = all the free electrons of a uniform planar foil have 
one and the same initial velocity Vq perpendicular to the 
foil. Rigorous results are obtained for test ion accelera- 
tion in the outer laminar zone (for strict definition see 
section ITlI Al below) of the dynamically evolving electron 
sheath. Particular attention is paid to the limiting en- 
ergy of accelerated ions at t — > oo, which is always finite 
within the adopted model. 

This paper is not the first publication addressing thus 
formulated problem: to a significant extent, it builds 
upon earlier work by Bulanov et al. [TtI ]. The new 
progress made here includes the following key issues. 
Temporal behavior of the boundary between the outer 
laminar and the inner relaxation zones of the dynami- 
cally evolving electron sheath is studied in detail. Con- 
sequently, the domain of applicability of the laminar-zone 
results for test ions is clearly identified in the full parame- 
ter space of the problem. In contrast to Ref. 01 > the elec- 
tric field in the laminar zone is calculated exactly and not 
to the accuracy of the linear in x term. It is proven that 
the answer to the intriguing question of whether a test 
ion can overtake the electron front (and, consequently, 
surpass the initial electron velocity vq) is determined by 
a critical relationship between the initial gamma-factor 
7o of accelerated electrons and the ion-electron charge- 
over-mass ratio [i [defined in Eq. below] . A fully rel- 
ativistic intermediate asymptotic formula (|58|l is derived 
which may be used in realistic situations to evaluate the 
maximum energy of accelerated protons and heavier ions. 



II. FORMULATION OF THE PROBLEM 

Consider a uniform electrically neutral plasma foil of 
thickness Iq with an initial density of free electrons no- 
At time t — all the free electrons are set in motion with 
the same initial velocity vq perpendicular to the foil (see 
Fig. [IJ; vo can be arbitrarily close to the speed of light 
c. At later times t > the motion of electrons, treated 
as a collisionless charged fluid, is governed by the electric 
field E(t, x), which arises due to charge separation in the 
evolving plasma cloud. The origin of the a;-axis, directed 
along the initial electron velocity vq, is chosen at the 
forward foil surface, so that initially the foil occupies the 
region — Iq < x < 0. The bulk foil ions are assumed to 
be infinitely heavy and staying at rest. Our goal is to 
calculate the motion of a test ion of charge +eZ p and 
mass m p placed initially at the foil surface x = 0, which 



is accelerated by the electric field E(t, x) in the positive 
direction of the x axis. Note that here and below m p is 
not necessarily the proton mass. 




FIG. 1: Schematic view of a plasma foil with motionless bulk 
ions (hatched area) and electrons (grey area) boosted to a 
velocity Vq. Dashed vertical line marks the initial position of 
electrons with the Lagrangian coordinate I. 

It is easy to understand that this problem is gov- 
erned by only three independent dimensionless param- 
eters, which we choose to be 

m e Z p l u! Q 2\-i/2 rt\ 

M= A = — , 70 = (l-j3g) ' ; (1) 

here m e is the electron mass, 

=(^) (2) 

is the plasma frequency in the initial configuration, e is 
the positive elementary charge, and /3q — vq/c. If a pro- 
ton (or some other light ion) is chosen as a test particle, 
the charge-over- mass ratio /i is small, u < 1/1836 <C 1. 
We will, however, explore the entire possible range of 
< jU < oo, firstly, for the sake of completeness of the 
analysis, and, secondly, keeping in mind that light pos- 
itive particles — such as positrons, ir + or u + mesons 
— may in principle be created and accelerated under an 
intense laser irradiation. 

For a non-thermal electron cloud considered here the 
quantity voy/jo/uiQ plays a role of the Debye length. 
Then, the parameter A is the initial foil thickness in units 
of the Debye length. Of particular interest is the limit 
A = of a geometrically infinitely thin foil which, how- 
ever, has a finite number 

£ = n lo (3) 

of electrons per unit surface area. On the one hand, the 
effects of charge separation are the strongest in this limit 
for a given value of vq- On the other hand, it is in this 
limit that most of the exact results can be established 



3 



analytically; many of them are then straightforwardly 
extended to a more general case of A > 1. 

Parameter 70 is the relativistic gamma-factor of the 
accelerated electrons. In the non-relativistic limit of 
ft <C 1, when 70 ~ 1 + \0a, this parameter becomes 
irrelevant, and we are left with only two principal param- 
eters \i and A. On a par with (3q and 70, the dimension- 
less momentum ITo = /3o7o is used below as the principal 
kinematic characteristic of the accelerated electrons. 



III. MOTION OF ELECTRONS 

A. General notation and relationships 

Motion of electrons is described by a function x e (t,l) 
[or by a function x e (t, £)], where I (or £) is a Lagrangian 
coordinate in the electron fluid: x e (t, I) is the position 
at time t of an electron whose original position at t = 
was x e (0,0 = -I (the broken line in Fig. [TJ. The 
dimensionless Lagrangian coordinate < £ < 1 is defined 
as 



(4) 



Then, the front edge of the electron cloud is at £ = 0, its 
rear edge is at £ = 1. 

At any fixed time t one can invert x e (t, £) with respect 
to £ to obtain the inverse function £ = £(i, x) defined 
inside the expanding electron cloud. Without electron- 
electron collisions the function £(£, x) ceases to be single- 
valued with respect to x after some time — even if it 
were so initially. In this paper we use the term "laminar 
zone" for that region of the electron cloud where £(t,x) 
is single- valued (see Figs. |21 and below). In the re- 
maining "relaxation zone", where £(t, x) is multi- valued, 
electrons gradually relax to the equilibrium Boltzmann 
distribution. 

The function x e (i,£) is found by solving the following 
equations of electron motion 



1 dx„ 



IF 



c dt v/I+nf 

dtt e e 

—TT = E(t,X e ), 

at m e c 



(5a) 
(5b) 



where (3 e = v e /c, and IT = j3 e ^ e = f3 e (l - 01) 1/2 ; the 
Lagrangian time derivative dj 1 dt is calculated at a fixed £. 
As is well known, equations JSJ admit the energy integral 



m e c 



1 +00 
; Sy 7e (i,£)d£+i- J E 2 dx = m e c 2 Y< lQ . (6) 



The electric field E(t,x) is obtained by solving the 
Poisson equation, which in the planar geometry of our 
problem yields 



here < a e (t,x) < 1 [0 < <Ji(x) < 1] is the fraction of 
the total number of electrons (ions) above x at time t, 
i.e. 



1 



a e (t,x) = — / n e (t,x')dx' 



(8) 



For motionless ions one obviously has 

1, x < -l , 
Oi(x) = { -x/Iq, -la < x < 0, 
0, x > 0. 



(9) 



Calculation of a e (t,x) depends on whether x happens to 
be in the laminar or relaxation zone of the electron cloud. 
In the laminar zone one simply has o e (t, x) = £(£, x), and 
equations of motion |[5j) can be solved analytically. 

In the relaxation zone, where the function £(i, x) be- 
comes multi- valued with respect to x, one has to use a 
more general expression 



a e {t,x) 



(10) 



where summation is done over all the segments A£ which 
at time t are located inside the interval (x, +00). Clearly, 
if Eqs. (0 are to be solved with a full account of the re- 
laxation zone, it can only be done numerically. To this 
end, a numerical code TIAC (test ion acceleration) has 
been written, which calculates the electron motion and 
the ensuing electric field E(t, x) with a full account of 
possible mutual interpenetration of different elements of 
the electron fluid. Because of rapid randomization (for 
details see subsection IIII B 2() of the electron motion at 
the core of the relaxation zone, such straightforward cal- 
culations can only be realized within a limited time span 
of the order of 100 periods of electron oscillations near 
the foil surface. 



B. Solution for A = 

Here we consider the limit Zo — > of a geometrically 
very thin foil with a finite value of electron number per 
unit area S = nolo, where both the initial electron den- 
sity no and the plasma frequency uo become formally 
infinite. In this case it is convenient to introduce the 
following units of time and length 



[x] = l s = v Q t s , (11) 



E(t, x) = 47reE [a e (t, x) — <Ti(x)] ; 



(7) 



[t]=ts = WE 70 ' 



which replace the usual time and length scales ^/to/woi 
voy/jo/ujQ in a finite-density plasma. Evidently, the 
length unit Iy, plays a role of the Debye length in our dy- 
namic electron sheath. Below, the quantities measured in 
units 1)1 1H are marked with a bar. The electron trajecto- 
ries are represented by a two-parameter family of curves 
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FIG. 2: Electron trajectories for A = in the non-relativistic 
limit of IIo <C 1. Each curve is marked by the corresponding 
£ value. The relaxation zone is shown as a grey shaded area. 

In the laminar zone, where <j e — £, Eqs. (J5J are solved 
analytically. The two relevant branches of this solu- 
tion, obtained with the initial conditions 5 e (0, £, /3q) = 0, 
n e (0,C,/3o) = n = /?o7o, are 



n e = n (i-£t), 
_ = t(2 - ff) 
" i + 

for < i < 2/£, and 

n e = no[(i-Ot + i-2/e], 



P (l-Q + 2t(l -2/£) + 4/£ 2 



(12a) 
(12b) 

(13a) 
(13b) 



for 2/f<t<2/K(l-0]- 

In the upper half-space x e > Eq. I|12b|) is easily in- 
verted with respect to £, which leads us to the following 
expression for the electric field 



E(t, x) = Aire's x) = 87re£ 



ftx 



(14) 



Inside the relaxation zone, where the electron trajec- 
tories intersect with one another, one has to abandon 
Eqs. (|12Jl - l|14|l and solve Eqs. (jHJ) numerically to calcu- 
late x e (i,£,Po) and E(i,x). 

In the non-relativistic limit x e (t, £) = x e (i,£,0) is a 
universal function of two variables i and £. It is plotted 
in Fig. [21 for a selection of £ values as calculated with the 
TIAC code. In Fig.[3]the non-relativistic function x e (i, £) 



is plotted versus £ for t = 20. Qualitatively, the relativis- 
tic trajectories x e (f, £,/?o) look similar to x e (i, £,0). In 
particular, for any < fln < 1 the rear edge of the elec- 
tron cloud £ = 1 turns around (i.e. has II e = 0) at i = 1, 
and later crosses the foil at i — 2 with II e = — 1. During 
the period < t < 2 there exists a vacuum gap between 
the foil ions and the ejected electrons. After its closure, 
the front and the rear edges of the electron cloud propa- 
gate freely in opposite directions with constant velocities 
±v . 
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FIG. 3: Position of individual electrons x e (t, f) in the non- 
relativistic case at time i — 20. Without electron-electron 
collisions, x e (i, £) is non-monotonic with respect to the La- 
grangian coordinate £ in the relaxation zone (shaded area). 
The number of local minima and maxima along £ rapidly (ex- 
ponentially) increases with time. 

Numerical solution of Eqs. JSJ enables one to trace the 
onset and subsequent expansion of the relaxation zone, 
shown in Figs. [21 and |31 as shaded areas. In this zone os- 
cillations of electrons around the positively charged foil, 
occurring on a time scale is, rapidly become stochas- 
tic and ultimately lead to establishment of the Maxwell- 
Boltzmann distribution with a certain temperature that 
can be found from the energy integral The starting 
point of the relaxation zone can be calculated analytically 
by applying the conditions dx e /d^ = 0, d 2 x e /d^ 2 = to 
Eq. (|13b[) : in the non-relativistic limit it has the coordi- 
nates 



t r o 



27 

T' 



r0 



which become 



4 
9' 

1 

2' 



X r (j 



X r0 



27 

32' 



(15) 



(16) 



in the ultra-relativistic limit 70 1. 



2. Evolution of the relaxation zone 

Even though electrons are treated as a collisionless 
fluid, the presence of a positively charged ion sheet is suf- 
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ficient for the electron subsystem to develop a stochasti- 
cal behavior. Once the electron trajectories in the central 
region begin to intersect one another in the process of os- 
cillations across the ion layer (see Fig. [21 , their further 
motion becomes increasingly stochastic. As a result, an 
isothermal quasi-equilibrium core gradually develops in- 
side the relaxation zone, where the electron density obeys 
the Boltzmann relation n e {x) cx exp (e(j>/T); here 4>(x) is 
the equilibrium electrostatic potential, and T is the tem- 
perature of the quasi-equilibrium core. This qualitative 
picture is fully confirmed by the results of direct simula- 
tions with the TIAC code presented in Fig. 0] 



TIAC: t= 100 
Boltzmann 




2 3 
ln(1 + 0.75 x) 

FIG. 4: Spatial profile of the electron column density a e (i,x) 
as calculated numerically with the TIAC code in the non- 
relativistic limit for i — 100. The equilibrium profile )[i 
obtained by using the Boltzmann relation is shown as a thick 
grey straight line. 



calculated from Eq. I|17|l . with that obtained from the 
TIAC numerical simulations in the non-relativistic limit 
for t — 100. One clearly distinguishes a quasi-equilibrium 
core of the relaxation zone, which is adequately described 
by the Boltzmann relation. Departures from the equilib- 
rium are significant in the outer part of the relaxation 
zone and, of course, in the laminar zone. 

To assess practical applicability of our results for test 
ion acceleration in the laminar zone, obtained below, one 
needs to know how the boundary x r (t) between the relax- 
ation and the laminar zones evolves in time. We analyze 
this evolution by combining an analytical estimate in the 
limit of t — > oo with the TIAC simulations for times 
t< 100. 

Since randomization of the electron motion in the re- 
laxation zone leads to establishment of the Maxwell- 
Boltzmann distribution, we can invoke the following ar- 
gument to evaluate the width x r of this zone at times 
( » 1: by a time i the relaxation zone spreads to a 
distance x r such that t = N r i eq {x r ), where N r is a nu- 
merical factor of the order unity, and t eq (x r ) is the travel 
time between x = and x = x r in the equilibrium po- 
tential l|17l) of an electron whose kinetic energy vanishes 
at x — x r . In a sense, t eq (x r ) is a timescale on which 
the information about maxwellization of the velocity dis- 
tribution up to a certain limiting value v r is transferred 
from x — to the corresponding limiting distance x r in 
the infinite potential well 1)17(1 . This argument sounds 
perfectly reasonable for the quasi-equilibrium core, and 
may be surmised to apply to the entire relaxation zone 
as well. 

The time i eq (x r ) is given by an integral 



Having adopted the Boltzmann relation, one easily 
solves the Poisson equation and obtains 



4>{x) 



-IT In 1 



m 

4T 



(17) 



where 4> and T are, respectively, the potential in units 
m e Ug7o/e and the temperature in units m e VQjQ. The 
temperature T is found from the energy integral JfjJ, 
which in our case transforms to a transcendental equation 



2Tpl + 



K (i/PqUqT) 

7o#i (1/ATIoT) 



= i; 



(18) 



here Kq(z) and Ki(z) are the modified Bessel functions. 
The limiting values of T are 



T 



n « i, 



i 

3' 

~ n » i. 



(19) 



Figure 0] compares the equilibrium electron column den- 
sity [as defined by Eq. 



1 f x 

T e,eq[X)\ s>0 ~ 2 V + 4T 



(20) 



teq {Xr ) 



dx 
v e (x) ' 



(21) 



where the velocity v e (x) is found from the energy integral 



7e(s) - 1 - 2/5 n f In 



AT+X r 

Af + x 



(22) 



for an electron moving in the static potential (|17f) . After 
some algebra we obtain 



teq(x r ) 



AT 



y/ln(l+x r /4T) 



1 + 2(3 n fz 2 

yf\ + A>n Tz 2 



dz. 



In the asymptotic limit of x r 3> 



(23) 

1 Eq. (EH yields 



teq \Xr) 



OLpqX t 



where the numerical coefficient a eo is 



a weak function of IIq: a e 



Sir for flo <C 1, and 



= 1 for IIo ^> 1. This leads us to a conclusion that 
asymptotically the ratio x r (t) jt should approach a cer- 
tain constant value a r (or oscillate in a narrow range 
around this value), which is fully confirmed by the TIAC 
simulations for times i < 100. 
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time t 



100 



FIG. 5: Evolution of the relaxation zone (shaded areas) for 
A = and three different values of 70 in terms of the hyper- 
bolic variable \r [see Eq. normalized to 770 = arccosh7o 
as calculated numerically with the TIAC code. The upper 
and the lower branches of the Xr(t) curves correspond, re- 
spectively, to the boundaries in the upper (x > 0) and lower 
(x < 0) half-spaces. Dashed curve shows the test ion trajec- 
tory for -y = l and fj, = 0.075. 



Figure El shows the temporal dependence of the upper 
and lower boundaries of the relaxation zone in terms of a 
hyperbolic variable Xr — Xr (t) introduced via a relation- 
ship 



tanh Xr 



x r (t) 
ct 



-(*) 



t 



(24) 



The fraction of the electron cloud occupied by the up- 
per laminar zone is given by the difference 1 — Xr/i]o [ m 
the ultra-relativistic case the separation along the hy- 
perbolic variable x — arctanh(a;/ct) is physically more 
representative than the separation along x], where 770 = 
arccosh7o = arcsinhllo. If Xr(t) should reach the value 
7]q, it would mean that the relaxation zone has reached 
the electron front x — i and the laminar zone has van- 
ished. One clearly sees that typically the electron cloud 
is not dominated by the relaxation zone, whose frac- 
tion asymptotically approaches some constant value, and 
which has a tendency to shrink with the increasing 70. 
For 70 = 1 this fraction is a r = 0.25 ± 0.02. Therefore, 
it can be expected that there exists a sizable window 
in the parameter space where the test ion acceleration 
takes place either entirely or predominantly in the lami- 
nar zone. 



C. Solution for A > 

For a finite-thickness foil it is convenient to introduce 
a different pair of time and length units, 

[t] = At s = uj Q 1 ^Ffy, [x] = Al s = voojq 1 ^/^, (25) 

where the time unit [t] is based on the relativistic plasma 

— 1/2 

frequency ^o7o • I n these units the frequency and the 



amplitude of the electron plasma oscillations are of the 
order unity for any value of 70 > 1. Below, the quantities 
measured in units H25[) are marked with a tilde. Note 
that in these units the dimensionless thickness of the foil 
is Iq = A. 



1. Electron trajectories 

For subsequent analysis of the test ion motion we need 
the electron trajectories in the upper half-space x > 0. 
However, because for A > each such trajectory x e (t, I) 
starts inside the foil at —A < x < 0, we should solve 
Eqs. (J5J in this region as well, with the initial conditions 
x e (0,l) = —I, II e (0,f) = n . The required relativistic 
solution in the laminar zone at —A < x e < has the 
form 



x e (t,l,P ) = -I + h(t,fio), 



(26) 



where h(t, /3rj) is a periodic function of time t; for the first 
quarter-period it is implicitly given by the quadrature 



1 - ¥ie 



1-c 2 



|/?o 2 C 4 



:d( 



(27) 



Actually, this is a solution for a relativistic particle mov- 
ing in a quadratic oscillator potential, for which one has 
the following energy integral 



1 



v^+nf + -/3 n h 2 



7o- 



(28) 



From Eq. (|28|l one readily establishes that the full ampli- 
tude (in units of vquJq 1 \/7o) of plasma oscillations inside 
the foil is 2ho, where 



ho 



270 
7o + l 



1/2 



(29) 



In the non-relativistic limit, obtained by putting /3o = 
in Eq. (|27|l . we have h(t) = sini and ho = 1. 

Now, the trajectory of an electron I in the laminar zone 
at x > 0, which matches the solution (|26J) . I|27|l at a point 
t = t e ±, x e = 0, can be written as 



x e (t,l,P ) 

Q = Q(lM 

R = R{t,l(3 ) 

tel = t e i(l, Po) 



R 



P+JP 2 -P 2 IR 



n„ 

1 - 



Q-l(t-t el ) 



1 



2 72 



Pol 



l 2 + \p 2 o~l\ 



1 



(t-tel) \2Q-l(t-t el ) 



l M 2 



c 2 + l M 4 



(30) 

(31) 
(32) 

(33) 
(34) 

(35) 
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here the values of / are confined to the interval < I < 
min{A,/io}- The non-relativistic limit of this solution is 
easily recovered by putting /3 = 0. To obtain the electric 
field in the laminar zone 



E(i, x) — 47renoWow 1 'Yq l(t, x) 



(36) 



one needs the function l(t,x,[3o), which is the inverse of 
x e (t, I, Pq) with respect to I. When calculating the results 
discussed in section lTVBI such inversion of Eqs. i|30|l - l|35|) 
was performed numerically. 

Since the Lagrangian coordinate I belongs to the in- 
terval < I < A, a vacuum gap is formed between the 
ejected electrons and the ion layer in the case of A < ho 
for a limited time after the rear electron edge I = A 
reaches the foil surface at x = — similar to the case of 
A = (see Fig. |2J where such a gap exists at < t < 2. 
Throughout the gap, the electric field E is constant and 
/ = A. To avoid unnecessary technical complications, we 
exclude the interval < A < ho < \/2 from our consid- 
eration, i.e. we perform all the calculations cither for the 
limiting case of A = 0, or for A > h$ when no gap ap- 
pears between the ejected electrons and the foil ions. The 
reward is that we get rid of the dependence of / on the 
parameter A, i.e. /(£, x, (3q) is a smooth function of £ and 
x for all £ > and < x < £, and this function does not 
depend on A. In particular, this leads us to an important 
conclusion that the behavior of the electron trajectories 
(hence, the distribution of any other physical quantity) 
in the laminar zone of the electron cloud at i > does 
not depend on the parameter A for A > /iq. 



2. Evolution of the relaxation zone 

The general arguments on the evolution of the relax- 
ation zone formulated in section ITlI B 21 for A = remain 
valid at A > as well. In particular, these arguments 
lead us to the following estimate for the relative width of 
the relaxation zone 



Xr(t) 



Vo 



< a xr < 00, 



(37) 



where Xr(t) is defined in Eq. (|24[1 . and a xr does not de- 
pend on time £ but is a function of A and 70. Unfortu- 
nately, we have no rigorous proof that one can introduce 
an upper bound l|37|) valid at all times £ < 00. We have 
only been able to verify it numerically with the TIAC 
code within a limited range of A < 20, £ < 100. As an ex- 
ample, Fig. shows the evolution of the relaxation zone 
for A = 10 and 70 = 1, 10, and 100. 

The TIAC simulations clearly indicate that the rela- 
tive width of the relaxation zone a xr increases with the 
increasing A, and decreases with the increasing 70. In 
particular, in the non-relativistic limit of 70 — > 1 we have 
a xr > 1 for A > 5.45, i.e. for A > 5.45 the laminar 
zone completely vanishes after some time (see Fig. |SJ. 




40 60 
time t 



100 



FIG. 6: Same as Fig. |^|but for a finite foil thickness A = 10. 
Shown here is the relaxation zone in the upper half-space 
x > above the initial foil surface x = 0. Dashed curve is the 
test ion trajectory for 70 = 10 and /i = 0.025. 



However, it reappears already at moderately relativis- 
tic electron energies 70 > 4, and becomes quite broad at 
70 > 10-100, when we have a xr < 0.2-0.3. The latter 
implies that the laminar-zone regime of test ion accelera- 
tion should be particularly relevant to highly relativistic 
electron clouds with 70 > 100. 

In the previous subsection it was established that the 
electron trajectories in the laminar zone do not depend 
on A for A > ho. This, however, does not mean that 
the same should apply to the boundary between the 
relaxation and the laminar zones because this bound- 
ary is determined by electrons traversing the relaxation 
zone. Nevertheless, one might expect that the curve Xr(£) 
should approach a certain limiting form for a fixed 70 and 
A — > 00: such a limit would correspond to the case where 
the bulk plasma ions occupy the entire half-space x < 0. 
However, the existence of such limit for all £ < 00 is 
far from obvious. Numerical simulations show that only 
the initial portion of the \r (£) curve around its first lo- 
cal maximum at £ ~ 15-20 (see Figs. |5] and 0) becomes 
independent of A for A > 2.5. 



IV. MOTION OF TEST IONS 



Acceleration of a positive test ion with a charge +eZ p 
and a mass m p (in conventional units) is described by 
the following equations of motion 



1 dx p 
c dt 

dn p 

~dT 



n 



V — 



n 2 

p 



V 



eZ r , 



E(t,x p ), 



(38a) 



(38b) 



where x p (t) is the position of the ion, and n p (£) = (3 p "f p = 
Pp(l— /3p) -1 ^ 2 is its momentum in units of m p c. As a rule, 
we assume that the test ion starts at £ = £ p o = with 



8 



the initial values x p (t p o) = H p (t p o) = 0. Then, because 
the electric field E(t,x) is non-negative for all x > 0, we 
are guaranteed that x p (t) > for all t > 0. Equations 
(|38|l bring in a dimensionless parameter /i = m e Zp/m p , 
which contains all the necessary information about the 
test ion. 



A. Solution for A = 

Here, as in section 1111 Bl we use the units 1)110 an d 
mark thus normalized quantities with a bar. The accel- 
erating electric field E(i, x) in the laminar zone is given 
by Eq. 



1. Vacuum phase 

When A = 0, a test ion begins its motion by pass- 
ing through a vacuum gap, where it is accelerated by a 
constant field E = 47re£, until its trajectory 



Xp^vacii) 



1 



i + /i 2 ng<2 - 1 



(39) 



crosses the rear edge [£ = 1 in Eq. I)12b|) ] of the electron 
cloud at 



f = t 



P i 



2(1 + M7o) 
1 + 2^70 + fj, 2 



£pi Xp^ vac (tp\) (40) 



with a momentum II 



pi 



/illot 



pi- 



The latter values 



should be used as the initial conditions for further ac- 
celeration inside the electron cloud. 



2. Non-relativistic solution 

In the non-relativistic limit a general analytical solu- 
tion is easily found to Eqs. (|38|) . The type of this solution 
depends on whether the roots A and 1 — A of the charac- 
teristic equation 



A 2 - A + 2/i = 



(41) 



are real or complex. In this way a critical value /i* = g 
of the parameter /i is established, which separates the 
two solution types. For /i > |, when A is complex, a 
test ion always catches up with the electron front x = t 
and acquires the final velocity v poo in excess of the initial 
electron velocity vq. 

In the physically more interesting case of /x < g (a suffi- 
ciently heavy test ion) the sought for solution of Eqs. 
is given by 



C 2 t\ 



where 



i(i-yr^) 



(42) 



(43) 



is the smaller of the two real roots of Eq. i)41)l . and the 
integration constants 



Ci 



2 A (1- X-fi) 
(1-2A)(1 + M )i+*' 



,1-A 



c 2 



(A-M) 



(l-2A)(l + /i) 2 - A 
(44) 

are calculated by using the initial conditions 1)40(1 for 70 = 
1; Eq. O applies at t > t pl = 2/(1 + ^). 

Solution ()42)) reveals the following general features 
of ion acceleration in the laminar zone. The distance 
i — x p (t) between the electron front and the acceler- 
ated ion increases monotonically with time, i.e. a test 
ion with H < \ lags further and further behind the elec- 
tron front as t — * 00. At the same time, the ion velocity 
v p {t) = dxp/dt monotonically grows in time, and in the 
formal limit of t — > 00 it asymptotically approaches the 
initial electron velocity vq — as it has been established 
earlier in Ref. [l7||. The latter, however, occurs on an 
extremely long timescale t ac ~ exp which is be- 

yond any realistic value for protons and other ions with 
M ^ 1/1836. From practical point of view, an intermedi- 
ate asymptotics 



, .'27re 2 E 
Vp(t) « 2[xvq In I 1 



m e v 



(45) 



inferred from Eqs. (02l"El|) for ^ < 1 and 1 <C i <C 
exp(i/i _1 ), might be of interest — if not the presence 
of the relaxation zone. The fact is that solution 142)1 ap- 
plies only for ions with 0.0745 < /i < 0.125 because, as 
one finds from the TIAC simulations (see Fig. [SJ), the 
ion trajectories for A = and /i < 0.0745 penetrate into 
the relaxation zone. Ions with fj, < 1/1836 are acceler- 
ated deep in the relaxation zone, where one can use the 
Boltzmann relation and the quasi-static potential 117)1 : 
one readily verifies that the potential ()17H leads to much 
higher (by roughly a factor /i -1 ) final ion energies than 
those obtained from Eq. ()45)l . From this we conclude that 
the acceleration in the near-front laminar zone of a non- 
relativistic electron cloud is never important for protons 
and heavier test ions. In section IIV Bl wc demonstrate 
that this conclusion, proved here for A = 0, is valid for 
A > 1 as well. 



3. Relativistic solution 

Analysis of ion motion in the general relativistic case 
is significantly simplified after we make a transformation 
from the dynamic variables x p and n p to hyperbolic vari- 
ables x an d f] defined by means of the relationships 

x x 

tanhx = — = Po^S-, sinh?y = U p . (46) 

Accordingly, the initial electron velocity Vq is represented 
by a parameter 770, where IIo = sinhryo, 70 = cosh 7/0, 
Po = tanh7/o. In terms of these variables equations 1)38)) 
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with the expression (|14|) for the electric field become 



dx 
dr\ 



coshx . , , 
cosh^ Smh(7? - X) 

„ coshx . . / x 

2/i — smh(?7o - x), 

cosh 77 



(47a) 
(47b) 



where £ = Int. Equations l|47|l do not contain the inde- 
pendent variable £ on their right-hand sides, which means 
that certain key features of the ion motion can be ana- 
lyzed by inspecting the integral curves of the first-order 
phase equation 



dq = 2 sinh(r/ - x) 
dx sinh(?y - x) 



(48) 



in the (x, v) plane. Of principal importance here is the 
singular point (x,v) = (Vo,Vo)- 

First of all note that physically meaningful in our con- 
text are the integral curves of Eq. I|48|) that lie in the 
half- plane 77 > x : this follows from inequality x p (t) < v p t 
valid for any motion with a; p (0) > and a monotonically 
increasing velocity v p (t). The electron front x = v$t is 
represented by the vertical line x = Vo- A test ion can 
cross the electron front either in a regular way at 77 > t/q , 
or by passing through the singular point (770,770) along 
one of the two characteristic directions defined by the 
characteristic equation l|41|) (see Fig. [7J. A remarkable 
fact is that the roots of Eq. I|41|) do not depend on 770, 
i.e. are the same for the relativistic and non-relativistic 
motions. As a consequence, we obtain a universal critical 
value //* = I which separates two topologically different 
patterns of the ion trajectories near the singular point 
(Vo,Vo)- 

For light ions with \i > i, when the singular point 
(%!%) is a focus (see Fig. |T^i), the qualitative picture 
is the same as in the non-relativistic case: a test ion 
always reaches the electron front x = Vo within a finite 
time interval and crosses it at 77 > 770 , i.e. with a velocity 
v p = v poo > v . 

For heavier ions with /i < | the singular point (x, v) — 
(770,770) is a node with two entrance directions 

77 — 770 = A(x — ?7o)j separatrix S, (49a) 
V ~ Vo — (1 — A)(x — Vo): general direction, (49b) 

where A is given by Eq. Note that for fi -C 1 we have 
A w 2/i <C 1. The separatrix S divides all the integral 
curves in the (x, rj) plane into two classes. To the first 
class belong the curves which lie below the separatrix S 
in Figs. and and enter the singular point along the 
general direction l|49b|l . From Eqs. igTJl, (gSJ and l|49b|) 
one readily verifies that these curves approach the elec- 
tron front x = Vo i n the asymptotic limit of t — > 00, with 
the value of 770 — x falling off as t _A . Exactly as in the 
non-relativistic limit, the distance t — x p (t) to the elec- 
tron front increases monotonically in direct proportion 
to as i — > 00, while the ion momentum H p (t) mono- 
tonically grows in time and asymptotically approaches 





[i = 7/128 <[i~| 



+ -1 
-2 



FIG. 7: Phase trajectories for the ion equations of motion in 
the form (1471) in the vicinity of the singular point (x, rj) = 
(770,770) for two values of parameter 77. The singular point 
is a focus for /i > 77, = | (a), and a node for fi < /i» (b). 
The separatrix S is a solitary integral curve of Eq. 14811 which 
enters the node (770,770) along the direction 77 — 770 = A(x — Vo), 
where A is given by Eq. 11431 . 



flo on a timescale t ac ~ exp (A -1 ). To the second class 
belong the trajectories that lie above the separatrix S 
and cross the electron front X — Vo at 77 > 770. Ions mov- 
ing along such trajectories overtake the electron front 
within a finite time and reach the final velocity v poo > v 
(ipoo > To)- 

A qualitative difference between the relativistic and 
non-relativistic cases arises when one considers the be- 
havior of the separatrix S. In the non-relativistic limit 
the latter is a straight line given by Eq. (|49a|l . There- 
fore, all physically interesting ion trajectories that start 
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X-Tl 




'crO 



A = 0, (j. = 7/128 
■q „ =9.0164 

'crOO 

Tl =14.131 

•crO 



-10 



FIG. 8: Global view of the relativistic phase trajectories of a 
test ion with fj, = 7/128. Depending on the ion start delay 
tpo > 0, the starting point of each trajectory may lie any- 
where between its intersection with the grey thick curve O, 
defined by Eqs. (I50L and with the bisector rj = \ (dashed 
line). Straight dotted segments of the phase trajectories cor- 
respond to ion acceleration inside the vacuum gap. 



from x p = with the zero initial velocity v p — at any 
time t p o > always lie below S, i.e. belong to the first 
class described in the previous paragraph. The relativis- 
tic separatrix, in contrast, bends down and crosses the 
line )j = x at a certain value r\ — rjo = X ~ Vo = ~Vcr00 
(see Pig. El, where ?7 cr oo = Vcrooin) is a function of /j. 
As a consequence, for a given /j and 770 > ?7crOo(/t) the 
phase trajectory of a test ion may pass above the sep- 
aratrix and fall into the second class. In such a case a. 
test ion overtakes the electron front within a finite time 
interval. To reach more definite conclusions, we have to 
take a closer look at the initial conditions. 

If a test ion begins to move simultaneously with elec- 
trons at t = t p o = 0, the initial part of its trajectory lies 
in vacuum and is represented by a segment of a straight 
line 77 = r/ vac (x) = 2% [as it follows from Eqs. I|39[) and 
|@BJ] with < % < Xpi — 5%>i; m Fig- El these segments 
are shown as dotted straight intervals. Then, the initial 
conditions for the phase equation (|48|l are given by the 
values 



1 

X = X P i = ^^pi' 



V = Vpi 



arcsinh 



2^sinhr/o (1 + /1 cosh rjo) 
1 + 2/i cosh r/o + fi 2 



(50a) 
(50b) 



inferred from Eqs. igOJ and igBJ. If we fix /j and 

treat 770 as a free parameter, Eqs. (|50() define a univer- 
sal curve O, the locus of the initial points for the integral 
curves of Eq. in the (x~ rjO: V~ Vo) plane (see Fig- - 
The value of parameter 770 along the curve O at its inter- 



section with the separatrix S defines the primary critical 
value r) cr Q = 77 cr o(/i) for this parameter. Its meaning is as 
follows: for any 770 > 7y cr o(/i) a test ion with the charge- 
over-mass ratio [i finally catches up with the electron 
front and overtakes it. 

The fact that the vacuum segments of the phase trajec- 
tories in Fig. El rise less steeply than the initial portions 
of the relativistic integral curves of Eq. (|48|) implies that 
a delayed (at t p o > 0) start of a test ion may result in 
its more efficient acceleration. In reality such a delayed 
start may occur when a test positive particle is created 
on the spot some time after the laser pulse. If we con- 
sider a delayed start at i p o > 2, when the vacuum gap 
is already closed (the opposite extreme to the previously 
considered case of simultaneous start at t p o ~ 0), the 
locus of the initial points (x, v) — (0, 0) for the integral 
curves of Eq. I|48|l in the (x — 770, V — Vo) plane will be the 
bisector line 77 — 770 = x ~ Vo- Hence, the intersection of 
the separatrix S with this bisector defines the secondary 
critical value 77 cr oo = ^ctOo(m) < Vcro(fJ-) of parameter 770 
(see Fig. EJ which has the following meaning: for any 
Vo < Vcroo (a 4 ) a test ion with the charge-over-mass ratio 
/_t always stays behind the electron front for all possi- 
ble starting times t P Q > . In the intermediate case 
of ?7 cr oo < ^0 < Vera a test ion can either overtake the 
electron front or stay behind it, depending on the start 
delay < i p o < 2. A selection of rj cr o(/j,) and 77 cr oo (/■*■) 
values calculated by solving Eq. I|48() numerically is given 
in Table |U 



TABLE I: Critical parameters rj cr o(n), T7 cr oo(/^), an d VcrA{n) 
as calculated for a selection of n values by numerically inte- 
grating the ion equations of motion 



1//1 


VcrO 


77crOO 


TjcrA 


10 


5.94216 


4.51193 


5.35497 


20 


15.7848 


9.88926 


15.0917 


100 


94.3602 


49.9827 


93.6670 


200 


193.685 


99.9916 


192.990 


1000 


992.089 


499.998 


991.395 


2000 


1991.40 


999.999 


1990.71 



4. The ultra-relativistic limit 

In the physically important limit of fj, <C 1 the func- 
tions r] cr o(ii) and r] cr oo(jj,) can be calculated analytically 
by using the ultra-relativistic (/1II0 3> 1) expansion of 
Eqs. (|5(3) for the curve O in Fig. El 



77 -770= In /i, (51) 

and the integral 

exp(77-?7o)-2 ) u(x-77o) =C (52) 

of Eq. in the limit of 770 - x > 1> V ~ X > 1 , 
when the right-hand side of Eq. I|48|l can be approxi- 
mated as 2/iexp(77o — 77). Having set the integration con- 
stant C = 1, we obtain the equation of the separatrix S 
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in the (x — fyo, ?7 — Vo) plane. Note that, although derived 
in the limit of rjo — X ^* 1) this equation has a correct lim- 
iting behavior at \ — > r/o as well. After we calculate the 
intersection points of the separatrix S with the curve O 
[as given by Eq. lf5T)l] and with the bisector x — V> we 
find 

?7 cr o = ln/i + — - 1, 7 cr0 = t:M ex P ( — - 1 ) ,(53) 
fi 2 VM / 



1 



?7cr00 = 7cr00 = ^ ex P . „ 

2/x 2 \2/j 



1 



(54) 



The corresponding critical values jcroin) and 7crOo(A t ) 
of the parameter 70 are obtained by applying the ultra- 
relativistic formula 7 = cosh 77 w | exp 77. Comparison 
with the numerical results from Table [I] shows that for 
> 20 the asymptotic formulae JHSJ), l|5*4*| for r] cr0 and 
Tycr-oo are accurate to within 1.4%. 

Making use of the integral l|52|l with the values of C = 
exp (rjoo — r/o ) > 1, we calculate the limiting value j poc = 
2 exp ijoo of the ion gamma-factor in the case when the 
ion overtakes the electron front, 



M7o 



l + ln 



270 



t p0 = and 70 > j C ro(fj), 
2/J7o ln(27o), i p0 > 2 and 70 > 7 cr oo (/•«)• 

This our result for 7 poo differs significantly from the value 
7 P oo = 27g calculated earlier in Eq. (35) of Ref. [l7| . 
which we believe to be erroneous. It should be noted, 
however, that Eq. 1155[1 can hardly be of any practical 
interest for protons and heavier ions because for /i < 
1/1836 the corresponding values of j cr o > 2.3 x 10 793 
and 7 cr oo > 2.4 x 10 398 are way too large to be ever 
encountered in nature. 



5. Intermediate asymptotics for the ion energy 

Having established that in reality test ions of common 
interest, i.e. those with /i < 1/1836, always stay behind 
the electron front, and that their dimensionless momen- 
tum H p (i) approaches the electron value IIo on an un- 
realistically long timescale t ac ~ exp (^/J -1 ), a natural 
step would be to look for an intermediate asymptotics 
for H p (t), valid at 1 <C t -C exp (^A* -1 ), that might be of 
practical interest for the problem considered. 

Once we let /1 -C 1 and agree that 70 <C 7ctOo(m) < 
7cr0 (/•*)% we can integrate Eq. (|47b|l in the limit of t ^> 1 
by making an approximation r\ — x, which is valid to the 
first order in /i along the direction of general approach 
(I49bll to the singular point (x,v) = (voiVo)- With the 
initial condition 

»7(Cpi) = Vpi, sinhr/p! = Il pl = (J,U i p i, (56) 

where £ p i = In t p \ and t v \ is given by Eq. (|4(J|) , the result 
of this integration reads 



tanh 



tanh 



exp[-2 M (C-C P i)]- (57) 



Performing Taylor expansion of Eq. I|57() with respect to 
the small parameter < 2/j(£ — Cpi) "C 1, we derive the 
following asymptotic expression for the test ion momen- 
tum H p {t) 



7 P 



= 1 



7o 



7 P i 



n 



pi 



7 P i + iipi \7pi + ripi 70 + n ) r ~ ' t p i 



7pi 

1/2 



fi In — 



(58) 

where 7 P = (11^ + l) i; ' and j pl = (ITp^ + l) 1/2 . Note 

that, once 1 < [< exp i^ 1 ), Eq. I|58|l applies at any 
degree of relativism of either electrons or a test ion, i.e. 
any of the three quantities IIo , n p i , and n p is allowed to 
be arbitrarily small or large compared to unity. Compar- 
ison with numerical integration of Eqs. I|47|l shows that 
for fi < 1/1836 the error of the intermediate asymptotics 
(JSBJ at 1 < ln(i/<pi) < 50 is typically about 2-4%, and 
never exceeds 8%. In both the non-relativistic (IIo *C 1) 
and the ultra-relativistic (/illo ^> 1) limits Eq. (|58|l re- 
duces to a simple expression 



n„ re n, 



pi 



1 



ln-t 



(59) 



B. Solution for A > ho 

From practical point of view the case of A > 1 is gener- 
ally more important than the limit A = considered so 
far. As a typical example, it may be noted that a 1 /im 
foil of solid gold ionized to z — 50 with 70 = 100 would 
correspond to Aw 10. Here we prove that practically all 
the qualitative and many of the quantitative results ob- 
tained for A = extend to the case of A > as well. To 
avoid unnecessary mathematical complications, we ex- 
clude from our consideration the intermediate range of 

1/2 

< A < ho and assume that A > ho = [270/(1 + 70)] • 
In units l|25|l the equations of ion motion H38fl become 



dx p 



dn 



EL 



1 + II2 



— ~- = fMal(t,x p ,po)- 
at 



(60a) 
(60b) 



As discussed in section IlII CI the function l(t,x,f3o) does 
not depend on A at A > ho, and is a smooth function of 
its arguments in the entire region 0<i<cx),0<i;<i. 
In the case of a simultaneous ion start at 7; p o = we have 
the initial conditions i p (0) = Hp(0) = 0, and the leading 
terms in expansion of the desired solution to Eqs. H60|) 
near t = are given by 



x p (t) 
IL p (t) 



^Wo? + O(t 4 ) 



(61a) 
(61b) 



Since neither Eqs. 1)60(1 nor the pertinent boundary con- 
ditions depend on A, we arrive at an important conclu- 
sion that, when expressed in units l|25|l . all the results 
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concerning test ion acceleration in the laminar zone are 
independent of A for A > h . 

The next important point is that in the limit of t — > oo 
equations (|fc>OI> become exactly equivalent to the corre- 
sponding equations of motion for x p (t), Tl p (t) in the case 
of A = 0. To verify this, we note that Eqs. [|3T))l - [f3T)|) 
imply I ~ 2/t <C 1 for f » 1. Then, by expanding 
Eqs. (|3U[) - I|35|) in powers of I, we derive an explicit for- 
mula 

l(i,x,Po) « 4(t-x) |P - Plx 2 + 2 + 

[(P $5? - 2) 2 + 8(1 - Pl)ix] 1/2 y] (62) 



which is valid for any < /3o < 1 in the limit of . 
and which is further simplified to 



2(t- x) 
t 2 - ffix 2 



« 1, 



(63) 



for t ^> 70 and all < x < t. Substituting Eq. I|63|) into 
Eq. (|60bjl . we obtain exactly the same equations with re- 
spect to x p (t), H p (t) as those with respect to x p (t), U p (t) 
in section llV Al which are then reduced to the same phase 
equation (|48fl with the same topology of integral curves 
in the vicinity of the singular point (XiV) = (VOtVo)- As 
a consequence, we arrive at conclusions that are fully 
analogous to those made for the case of A = 0: 

(i) test ions with /i > /i» = g always catch up and over- 
take the electron front, being finally accelerated to 
n poo > IIo; 

(ii) for any fj, < g there exists a critical value ?7cr-A(/-0 
of parameter 770 [or, equivalently, a critical value 
7crA(/-0 °f parameter 70] such that only for 770 > 
ricrh(n) [i-e. for 70 > TwaO-*)] can a test ion with the 
charge-over-mass ratio fi overtake the electron front 
and be accelerated to TL poo > IIo; f° r ^0 < VcrAifj) 
a test ion lags behind the electron front while its 
dimensionless momentum lip (t) asymptotically ap- 
proaches n on a timescale t ac ~ exp (^M 1 )- 

These conclusions are fully confirmed by numerical in- 
tegration of Eqs. H6Uf) with l(t,x,/3o) calculated from 
Eqs. 

A remarkable fact is that we have two universal func- 
tions, namely, rj cr o(fj,) for A = 0, and rj cr \(fj,) for A > ho, 
which cover the entire range of A variation. Since at 
intermediate times t ~ 1 the two cases of A = and 
A > ho are mathematically not equivalent, the functions 
Vcro(lJ') and ijcrA^) numerically differ from one another, 
except for the initial point ^ cr o(M*) = "Her A (M*) = 0. This 
difference, however, is practically not significant, as one 
verifies by comparing the numerically calculated values 
of rj cr Q(^L) and rj cr \{i£) in Table [I] In the limit of /i < 1 
one derives an asymptotic expression 



1 ( l 1 

7crA — ex P 1 

4 



(64) 



which is exactly one half of the corresponding limit i|53|) 
for jcrQ. One can safely conjecture that for intermediate 
values < A < ho the corresponding critical values of 
70 lie between 7 cr0 (M) and 7 ci -a(m)- For 70 > 7crA(^) ; 
when a test ion does overtake the electron front, its final 
energy in the limit of fj, <C 1 is given by 



7poo = M7o 1 + In 



47o 



(65) 
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FIG. 9: Comparison of the intermediate asymptotics l|58|l 
with the results of numerical integration of the ion equations 
of motion 1601 for fi — 1/2000 and four different values of 70: 
plotted is the ratio of the ion momentum Y[ pap (t), obtained 
from Eqs. I58H and 166H . to the numerically calculated value 
Ilp(t) versus time t in units 11251 . 

Under the approximation l|63|) one derives the same 
intermediate asymptotics l|58|) for the test ion momentum 
Hp(£) as in the case of A = 0. The only difference between 
the two cases is in the initial values of t p \ and II„i — 



U p (t p i) to be used in Eq. (|58j|. Unlike in the A = 
case, no appropriate analytical solution to Eqs. (|60|l was 
found for t<l that would yield suitable expressions for t p \ 
and Upi. It is only in the ultra- relativistic limit /1II0 2> 
1 that one obtains a simple result Upi = /xllo, t p \ = 

(k/£7q) . In the opposite limit of /iHo <C 1 expansion 
(IGlbll suggests that the value H p ~ /iil is achieved at 
t ~ 1. Taking guidance from such considerations, we 
propose simple approximate expressions 



n p i = ^n , t pl = [ 2 



2^70/^0 



1/2 



(66) 



which, on the one hand, agree with the ultra-relativistic 
limit and, on the other hand, fit reasonably well the re- 
sults of numerical integration of Eqs. I|6U|) shown in Fig. [5] 
As a result, Eq. I|58|l with n p i and t p i taken from Eq. I|66|) 
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is formally applicable at t p i < ( < exp (i/i" 1 ). A com- 
parison with numerical results in Fig. [9] shows that at 
t > 100 the intermediate asymptotics Ij58(l has a typical 
error of a few percent. 

If we consider now the case of a delayed ion start at 
t = t p o > 0, we find that for t P Q > 1 we have I <C 1 for all 
t > t P Q, which again leads us to the approximation l|63() 
and to the phase equation (|48|) with the initial condition 
(x>v) — (0,0). Hence, exactly as in the A = case, 
the critical value of the 770 parameter for t p o 3> 1 should 
be given by the function rj cr Qo(^,). This conclusion is also 
fully confirmed by numerical integration of Eqs. I|60(l . For 
the dependence of the critical 770 on the ion start 
delay t p o is as follows: as t p o increases from t P Q = to 
t p0 = ^/2, the critical value of 770 decreases from rj cr \(fi) 
to ?y C r oo(M), an d for t p o > y/2 it remains equal to r; cr .oo(/x). 

C. Domain of applicability of the laminar-zone 
solution 

The foregoing analysis of test ion acceleration has been 
based on the assumption that the ion trajectories lie en- 
tirely inside the laminar zone of the electron sheath. As 
already mentioned in section lTTTl this is true only within a 
certain domain of our three-dimensional parameter space 
(/i,7o, A). Because of a weak dependence on A, the lim- 
its of this domain can be conveniently analyzed in the 
two-dimensional (/-t _1 ,?7o) plane. 

Let rji am (/x) be the threshold value of 770 = arccosh 70 
at which the trajectory of an ion with a given charge- 
over-mass ratio /i just touches the outer boundary of the 
relaxation zone (see Figs. [5] and EJ), i.e. for 770 > 77; aTO (/x) 
the ion trajectory lies entirely in the laminar zone, and 
for 770 < r)i am (/j,) it penetrates (at least partially) into the 
relaxation zone. Typically this touching occurs near the 
first local maximum of the corresponding x r (t) curve at 
t, t ~ 15-20 and lies within the reach of the TIAC code. 
Because this early part of the boundary between the two 
zones ceases to depend on A for A>2.5, the same applies 
to the function r]i am (n). 

Figure ^| shows two curves r)i am (fJ-), calculated for 
A = and 10, which actually span the entire dependence 
of 772am on the A parameter. The domain of applicability 
of the laminar-zone results lies outside the grey shaded 
area bounded by the 772am (m) curves. For protons with 
\i = 1/1836 it corresponds to 70 > 348 at A = 0, and to 
70 > 537 at A > 2.5. The fact that the two critical val- 
ues r) C ro(p,) and r) cr oo(fi) turn out to be deeply inside this 
domain (at least for /i < 0.1) justifies all the conclusions 
made in sections llV Al and IIV Bl about the possibility for 
a test ion to catch up with the electron front. Note that 
the two curves ri cr o(fi) and rj cr \(fi), which span the de- 
pendence of the critical value rj cr on the A parameter, 
are virtually indistinguishable in Fig. 1101 

For protons and heavier ions with fj, < 1/1836 the fol- 
lowing conclusions can be drawn from Fig. 1101 The posi- 
tion of the rjiam (a*) curves indicates that ion acceleration 




1 +ri 

'o 

FIG. 10: Characteristic regions in the (770, /i _1 ) parameter 
plane. Domain of applicability of the laminar-zone solution 
lies outside the grey shaded area delimited by the two 77z am (/i) 
curves, the inner one calculated for A = 0, and the outer one 
for A = 10. Hatched region is where test ions overtake the 
electron front: it is delimited by the ?7 cr o(/z) curve which is 
practically indistinguishable from r) cr ^{fi). Dashed horizontal 
line corresponds to the proton value fi -1 = 1836. 



in the non-relativistic case of Ho <C 1 takes place deeply 
in the relaxation zone, where one can expect the usual 
Boltzmann relation to be a good approximation. The 
laminar-zone solution is not applicable in such a case. 
However, it becomes fully applicable when the electrons 
are boosted to highly relativistic energies of 70 > 300-500 
fa) > 6-7). 



V. CONCLUSION 

In this paper rigorous results are presented for a partic- 
ular case of the TNSA mechanism of ion acceleration in 
a planar electron sheath evolving from an initially mono- 
energetic cloud of hot electrons. Self-consistent treat- 
ment of the collisionless electron dynamics fully captures 
the effects of departure from the Maxwell-Boltzmann dis- 
tribution. These effects come to a foreground in the outer 
laminar zone of the expanding electron cloud, where the 
ion acceleration can be analyzed by analytical means. In 
particular, the limiting (in the limit of t — > 00) gamma- 
factor 7 poo of an accelerated test ion can be calculated ex- 
actly. Note that the assumption of an isothermal Boltz- 
mann distribution for hot electrons leads to an infinite 
value of 7 poo . 

It is shown that the limiting value 7 poo is determined 
primarily by the values of the two (out of the total three) 
principal dimcnsionless parameters of the problem, the 
ion-electron charge-over-mass ratio ji = m e Z p /m p , and 
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the initial gamma-factor 70 of the accelerated electrons. 
For /x > /i* = i a test positive particle (for example a 
positron) always overtakes the electron front and reaches 
7poo > 7o- I n the physically more interesting case of 
fi < /1* the limiting ion energy depends on whether 70 
is above or below a certain critical value 7 cr = 7 cr (/z), 
namely, we have j poc = 70 for 70 < j cr , and 7 p0 o > 70 
[as given by Eqs. 1 (55)1 and J55J] for 70 > 7 cr . Practically 
insignificant dependence of 7 cr on the dimensionless foil 
thickness A is limited to a variation within a factor 2 and 
spanned by the functions 7 cr o(/i) and 7 cr A(^t) calculated 
in Table Hand Eqs. Q and (S3) . 

For protons and heavier ions with /i < 1/1836 we al- 
ways have 70 < 7cr because the corresponding values of 
7 cr ~ exp(/i~ 1 ) are enormous and beyond practical reach. 
Therefore, in reality these ions can never catch up with 
the electron front. Although formally the ion gamma- 
factor j p (t) in this case still tends to 70 as t — > 00, this 
fact is also practically irrelevant because "f p (t) approaches 
70 on an enormous timescale t ac ~ isexp(i/i _1 ) [or 
t ac ~ (jJq -1 ^/^) exp Qm - 1 ) f° r A > 1] that never occurs 
in nature. For practical applications one should use the 
intermediate asymptotic formula i|58|) derived for dimen- 
sionless times 1 <C i, t <§; exp (^M • 



Our results for ion motion have been obtained un- 
der the condition that the ion trajectory lies entirely in 
the laminar zone of the electron sheath. Numerical in- 
vestigation of the evolution of the laminar zone bound- 
aries reveals that this condition imposes a lower bound 
7o > Jiamip) on the initial gamma- factor of hot electrons 
(see Fig. I10J) . The latter inequality reflects a more gen- 
eral fact that the role of the non-Boltzmann effects in ion 
acceleration increases with 70, i.e. with the energy of hot 
electrons. In particular, acceleration of protons by a non- 
relativistic electron sheath occurs practically entirely in 
the quasi-Boltzmann core of the electron cloud, where the 
details of the initial electron energy distribution are "for- 
gotten" . However, in the ultra-relativistic case of mono- 
energetic electrons with 70 > 300-500, acceleration of pro- 
tons takes place entirely in the non-Boltzmann laminar 
zone of the electron sheath, where full memory of the 
initial electron energy distribution has been preserved. 
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